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ABSTRACT 

We herein determine the average integrated mass loss from stars belonging to the 
dominant metal-poor population ([Fe/H] ~ -1-7) of the Galactic globular cluster oj 
Centauri (NGC 5139) during their red giant and horizontal branch evolution. Masses 
are empirically calculated from spectroscopic measurements of surface gravity and 
photometric measurements of temperature and luminosity. Systematic uncertainties 
prevent an absolute measurement of masses at a phase of evolution. However, the 
relative masses of early asymptotic giant branch stars and central red giant branch 
stars can be measured, and used to derive the mass loss between these two phases. 
This can then be used as a physical check of models of horizontal branch (HB) stars. 
For oj Cen, the average difference is found to be 26 ± 4 %. Assuming initial and final 
masses of 0.83 and 0.53 M©, we determine that 0.21 ± 0.03 M Q is lost on the RGB 
and 0.09 ± ~0.05 M is lost on the AGB. The implied HB stellar mass of 0.62 ± 0.04 
M Q is commensurate with literature determinations of the masses of the cluster's HB 
stars. The accuracy of this measurement can be improved through better selection 
of stars and spectral coverage, and applied to other clusters where horizontal branch 
models do not currently agree. 

Key words: stars: mass-loss — circumstellar matter — infrared: stars — stars: 
winds, outflows — globular clusters: general — stars: AGB and post-AGB 
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1 INTRODUCTION 

Though all stars lose mass, the vast majority experience 
their most significant mass loss on the red and asymptotic 
giant branches (RGB/ AGB). The amount of mass lost de- 
termines both the mass and type of the stellar remnant. It 
can profoundly alter not only the post-RGB track of a star's 
evolution in a Hertzsprung-Russell diagram (HRD), but also 
the mass, chemical state, composition and mineralogy of the 
material returned to the interstellar medium. The processes 
that govern thi s mass loss a re relatively-well understood 
(e.g., see review IWillsonl lioOOj . and can be broadly divided 
into two categories. Mass loss from a star's hot (~10 000 
K) chromosphere eject s mostly atomic gas at low mass-loss 
rates (<1CT 7 M Q yr-^ lDupree et al.|[2009h over most of the 
star's life. Conversely, pulsation-enhanced mass loss driven 
by a cool (<2000 K), dusty wind can reach rates of >10 
Mq yr _1 (|van Loon et al.ll 19991 ). but for a much shorter pe- 
riod of time. 

Instantaneous mass-loss rates have been derived for a 
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large number of stars, but they are difficult to measure 
with any absolute accuracy and can vary substantially on 
cvolutionarily-short timescales. Determination of the inte- 
grated mass-loss rate (i.e. the difference in mass) between 
two evolutionary phases is necessary if one is to probe how 
a star with a given set of parameters will evolve. 

Stellar clusters provide excellent testbeds in which to 
probe stellar evolution, containing one (or sometimes a 
few) distinct populations at a known age, distance and 
metallicity, which have giant stars at the same initial 
mass. Increasingly-precise stellar ev olution models (e.g. 
iMarigo et alJl2q08l:lDotter et al.ll2008l ) and white dwarf ob- 
servations (e.g. Moehler et al. 20041 ) can determine the ini- 
tial and final masses of present giants to within a few per- 
cent. In intermediate stages, however, there are considerable 
systematic differences in the modelled stellar masses. 

This is perhaps most notable for horizontal branch (HB) 
stars. In clusters with blue HBs, stellar masses are relatively 
easy to determine via HB modelling and RR Lyrae pulsa- 
tions. If the masses of HB stars are >0.65 Mq (depending 
on the HB model and the cluster metallicity), they will be 
too cool to become RR Lyrae stars and modelling becomes 
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less precise due to the weaker dependence of the predicted 
Teff of red HB stars on their mass. In some cases, different 
models can disagree on HB ma sses by as much as ^25%, de- 
pending on the model chosen (ICarraro et al.|[l996l ; ICatelanl 
l2009bl ; [McDonald et al]|2011bl '). This causes severe problems 
when one tries to determine the total mass lost on the RGB, 
whether RGB or AGB mass loss is dominant, and the de- 
tailed evolutionary stages where mass loss is important. 

In this paper, we use a simple method to directly deter- 
mine the fractional mass difference between RGB and early- 
AGB stars, which by definition gives the time-integrated 
mass loss between these two phases. We apply this method 
to the metal-poor ([Fe/H] ~ -1-7) population of u) Centauri 
(NGC 5139): a well-studied globular cluster with a blue HB 
and well-determined HB star masses. This allows us to em- 
pirically confirm the accuracy of the HB models and inves- 
tigate the precision with which we can determine the time- 
integrated mass loss. 



2 METHOD 

If one knows the distance and reddening to a star, one can 
determine its mass simply from three observable param- 
eters: the gravity, log g, determined from high-resolution 
spectroscopy; the bolometric luminosity, L, derived from 
the integrated spectral energy distribution (SED) and clus- 
ter distance; and the effective surface temperature, T e g, de- 
rived from spectroscopy, photometric colours or SEDs. Us- 
ing these, the stellar radius can be simply derived using the 
Stefan-Boltzmann Law, and Newtonian gravity as follows: 



Table 1. Spectroscopic and photometric parameters for our tar- 
get stars. A complete table is available online. 
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Given the difficulty in determining these parameters accu- 
rately (particularly log g), it is not normally possible to use- 
fully constrain the masses of individual stars. Systematic er- 
rors mean that absolute masses for groups of stars are also 
usually too imprecise to be useful. By applying a statisti- 
cal approach to a large number of systematically-surveyed 
stars, however, one can gain a useful measure of the dif- 
ference in mass between two populations of stars within 
the same cluster that is theoretically model-independent. In 
practice, one must accept some dependency on the stellar 
atmosphere models used to derive the observables (L, T e g 
and log g), but these can mostly be circumvented by choos- 
ing stars with a similar temperature range where systematic 
effects are small. 

In recent years, large-scale studies of stars in globu- 
lar clusters have meant many stars have had their funda- 
mental parameters accurately determined. The cluster cj 
Cen has had photometric temperatures and luminosities, 
and spectroscopic metallicities determine d for a large num- 
ber of its stars dMcDonald et al.l 120091 . hereafter M+09; 
and [johnson & Pilac howskil |2010L hereafter JP10, re spec- 
tively: see also lvan Loon et al.ll2007l : lMarino et al.ll201ll ). We 
present here the spectroscopic gravity determinations from 
the original spectra of JP10 and use these to measure the 
relative masses of RGB and AGB stars. 
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3 SAMPLE SECTION & DATA REDUCTION 



An essential requirement is to accurately differentiate 
RGB from AGB stars. The most reliable way to do so 
is photometrically, on the early- AGB (eAGB), where the 
two branches separate on Hertzsprung-Russell and colour- 
magnitude diagrams (HRD, CMD). Our selection is based 
on cuts app l ied tp the HRD of M+09 and the CMD of 
iBellini et atl i|2009h : stars common to M+09 and JP10 are 
selected if they fall into the designated regions on both the 
HRD and CMD (Figure Q). 

A well-known spread exists in the elemental abundances 
and metallicity of a; Cen's stars (e.g. JP10). To minimise the 
effects of this spread, and due to the difficulty in identifying 
AGB stars in the metal-rich populations in the cluster, we 
further restrict our sample to stars with a metallicity range 
of -1.9 ^ [Fe/H] ^ —1.5. This is centered on the cluster's 
bulk population, at [Fe/H] = -1.7. This also helps prevent 
any metallicity-related bias, and leaves a total sample of 161 
RGB and 38 eAGB stars. 

Changing the metallicity limits we choose (-1.9 ^ 
[Fe/H] ^ -1.5) by ±0.1 dex does not change our results by 
more than the random error. We also probed for variations 
in the mass differential within the sample, but could find no 
significant variation with metallicity or with any other ele- 
mental abundance (as listed in JP10). Seven of the stars in 
our final sample (LEIDs 23033, 24027, 38226, 42174, 43104, 
43108 and 47151) have positive [Na/O] abundances suggest- 
ing they may belong to a second-generation, helium-rich 
population. These are all RGB stars, and excluding them 
makes no significant difference to the final result. 

A systematic temperature difference of ~130 K exists 
between the reddening-corrected photometric temperatures 
of M+09 and the (V-K )-based temperatures of JP10. This 
was traced to inaccuracies in the filter transmissions and 
zero point fluxes used for the optical photometry in M+09. 
The M+09 data were re- red uced, adoptin g a reddening value 
of E(B — V) = 0.12 mag (jHarria Il996l ) and a distance of 
5.3 kpc (see Sj4] below), an d replacing the ori g inal p hoto- 
graphic BV photometry of Ivan Leeuwen et al.l 12000) with 
the UBVRI CCD photometry of lBellini etHl (|2009l ) . This 
reduced the temperature difference between the two datasets 
to an average of 47 K (r.m.s. 56 K). We retain the sample 
selection used above as the M+09 and Bellini et al. data are 
independent samples, and the M+09 data shows less scatter 
on the HRD. 

The \ogg values published in JP10 are based on photo- 
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metric colours, which assume a priori a mass of 0.8 Mq. As 
we measure mass, we require spectroscopic values of logg, 
which are free from such assumptions. These were deter- 
mined by comparing the [Fe/H]i and [Fe/H]n abundance^ 
derived from the equivalent widths provided in JP10 (their 
Table 3). Model atmosphere grids were constructed for each 
star with logg covering ±0.5 dex (in 0.01 dex increments) 
from the photometrically-derived value provided in JP10 
(their Table 2). Although the T c g and microturbulence val- 
ues were held fixed at the values given in JP10, we iteratively 
adjusted the model metallicity to equal the average value of 
the [Fe/H]i and [Fe/H]n ratios for a given logg grid point. 
The final spectroscopic gravity was then selected as the log g 
value that satisfied [Fe/H]i = [Fe/H]n. The need for accu- 
rate logg values meant that spectra were discarded where 
one or both of the Fe II lines could not be reliably measured 
and/or the two lines did not provide reasonable agreement 
on the "best-fit" logg value. The comparative weakness of 
the Fe II lines (Figure [2]) means that this represents a sub- 
stantial reduction in our sample size. A further four RGB 
and two AGB stars were removed because they no longer fell 
into their respective regions in the HRD created using the 
Bellini et al. data. The final sample of high-quality targets 
comprises of 66 RGB and 21 eAGB stars. This subset of 
targets has negligible difference in their average metallicity, 
temperature and luminosity compared to our original selec- 
tion. The new photometric and spectroscopic parameters for 
these objects are listed in Table [1] where stars are listed by 
their Leiden Identifier (LEID). These parameters were then 
used to measure a mass for each star, using Eq. {TJ. 



4 ERROR BUDGET 

The random error associated with individual RGB stellar 
masses is dominated by errors in logg caused by uncer- 
tainties in the stellar metallicity, temperature and luminos- 
ity. The metallicity errors (taken from JP10) and a con- 
servative microturbulence error of ±0.3 km s _1 were prop- 
agated through the logg determinations described above. 
From these, a mass error is derived of SM/M = 1 10 %. In- 
ternal temperature errors of f»56 K were estimated (being 
the r.m.s. deviation between values in Table [T] and JP10), 
yielding an error of 5M/M = 1"]J%. Finally, luminosity er- 
rors of ~5% yield a proportional effect on SM/M. Taking 
these errors in quadrature, we derive a r.m.s. random er- 
+ 18 o/„ Errors on the AGB should be sim- 



ror of SM/M = ±\*% 
ilar, though the logg determination imparts a ~^\®% error 
in these stars, bringing the total error to SM/M = ±f %. 
This matches well with the r.m.s. scatter observed in Figure 
[3] the r.m.s. scatter is 20% for the RGB stars and 17% for 
AGB stars. This shows that our random error budget cov- 
ers the observed errors and suggests that the spread in the 
masses of stars on the eAGB is unlikely to be more than a 
few hundreds of a solar mass. 

The systematic error budget is dominated by errors in 
temperature, reddening and distance to cj Cen. Systematic 
errors in temperature of ~50 K were estimated, based on 



[Fe/H] abundances derived from Fe I and Fe II lines, respec- 
tively. 



600 
500 

400 

o 

m 300 
TS 

ii- 200 



100 



80 



V* 



5200 5100 5000 4900 4800 4700 4600 4500 4400 4300 
Temperature (K) 



12.6 
12.8 

13 
13.2 
13.4 h 
13.6 
13.8 

14 




0.7 



0.9 



1.1 



1.2 



1.3 



(B-V) 



Figure 1. Hertzsprung-Russell and colour- magnitude diagrams 
showing the central giant branches stars in ui Cen. The larger 
red and green circles show our RGB and eAGB targets, resp ec- 
tively. An isochrone and AGB model from Pott er et al.l i2008T ) is 
included, set at 12 Gyr, [Fe/H] = -1.7, [a/Fe] = +0.4. 



the 47 K systematic difference in the temperatures in Table 
[T]and those of JP10, providing an error of SM/M = ±11% 
as above. A conservative error of SE(B — V) = ±0.025 also 
yields a systematic temperature error of ~50 K (M±09), 
providing an identical error. Dista nce estimates to u ) Cen 
vary between 4.8 and 5.52 kpc (|van de Ven et al. I l200d : 
Idel Principe et all I2006T ). We therefore assume a distance 
error of Sd/d = tg% to encompass these values. This gives 
SL/L = SM/M = tig%- Combining the temperature, red- 



dening and distance errors in quadrature, we find a system- 
atic error of SM/M = ±H%. 

While difficult to quantify, additional sources of error 
may arise from the absolute accuracy of the MARCS model 
atmospheres used to measure the stellar luminosities and the 
wavelength dependence of the reddening correction. Choice 
of atmosphere models has very little effect: differences be- 
tween masses derived using the Kurucz, Kurucz a-enhanced 
and MARCS atmospheres are <1 %. Using (V — K) c olours 
to derive a bolometric luminosity (|Alonso et al.|[l999h . how- 
ever, yields masses which are ~5% greater than those de- 
rived from SED fitting. A significant difference here may 
be the effect of the H~ opacity minimum at 1.6 /im. The 
observed flux in the 2MASS H filter is ~6% less than ex- 
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Figure 2. Example spectra of an RGB and AGB star with a 
typical temperature and metallicity (LEID 56118, T c g = 4620 
K, [Fe/H] = -1.64; and LEID 34056, T cff = 4875 K, [Fe/H] = 
-1.63, respectively). Synthetic spectral models of the 6149 and 
6247 Fe II lines are overplotted with a fixed [Fe/H]ll abundance. 
Open circles indicate the observed spectrum, the solid black line 
shows the best-fit \ogg synthesis (derived from the EW analysis), 
and the dashed lines indicate changes in logg of +0.1/-0.1 dex 
(blue/red lines). 



pected (M+09). The accuracy of colour corrections to U 
and B band data in iBellini et all HH) may also be im- 
portant: Bellini et al. note that their U magnitudes may be 
systematically incorrect by up to 0.15 mag. Removing the 
[/-band data from our SED fits reduces the derived mass by 
~6%. Similarly, altering the reddening law within a reason- 
able range can produce errors of ~3%. Combining the above 
uncertainties yields an approximate error of SM/M « ±7%, 
giving a final systematic error of SM/M ~ %• 

Our sample stars cover a relatively-narrow tempera- 
ture regime, and the derived masses show no statistically- 
significant variation with temperature within the RGB pop- 
ulation (Figure[3|. The eAGB population shows a weak tem- 
perature trend, though this is based on small-number statis- 
tics may be intrinsic: stars with lower envelope masses will 
have hotter temperatures, as on the HB. While the afore- 
mentioned systematic errors will affect the derived masses 
of our sample, the relative masses of the RGB and AGB 
populations should remain largely unaffected. The ratio of 
the masses of the RGB and eAGB samples is therefore our 
most robust determination of mass loss between these two 
evolutionary points. 



5 THE RGB VS. EAGB MASS RATIO 

The average masses and their sample standard deviations 
found for our sample of RGB and eAGB stars, respectively, 
are 0.657 ± 0.090 M and 0.488 ± 0.117 M©. We stress that 
these values do not take into account the systematic errors of 
SM/M = ±2°% outlined above. The RGB mass is 18% lower 
than the c anonical ~0.8 Mq ex pected from stellar isochrone 
modelling (|Girardi et al.ll2002l : M+09), which is commensu- 
rate with these errors. It would be speculative to suggest 
why the masses we derive are systematically lower, though 
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Figure 3. Variation of derived mass with temperature for our 
RGB (hollow red circles) and eAGB (filled blue squares) samples. 



we note that the cluster's distance is (by a small margin) 
the largest contributor to the systematic error budget. 

Taking the ratio of the two masses, which is largely 
independent of systematic errors, we find that 25.7 ± 4.3 
% (standard error) of mass is lost between the central RGB 
and eAGB. The true error may be a little larger (perhaps 
1-2%), due to temperature-dependent effects below the level 
of statistical detectability, and the robustness with which we 
can photometrically differentiate RGB and AGB stars. 



6 DISCUSSION 

While the absolute masses we determine for our RGB and 
eAGB stars are systematically uncertain, the percentage dif- 
ference between them can be put into context using well- 
calibrated initial and final stellar masses. This allows the 
mass of the HB stars to be determined and compared with 
HB models. The mass predicted by isochrone fitting for a 
200 L RGB star i n uj Cen is between 0.80 and 0.85 M 
l|Girardi et all l2002i : iDotter et af] 120081 : M+09, and refer- 
ences therein), with very little mass loss occurring before this 
stage. This mass range agrees with the 0. 806 ± 0.056 M 
deter mined from the binary OGLEGC 17 l|Thompson et all 
l200ll) . The final white dwarf mass in similar globular clus- 
ters is found to be 0.53 ± 0.03 Mq (also derived using Eq. 
<rT)>: lMoehler et alj|2004h . 

Taking these two limits, this yields an eAGB star mass 
of 0.62 ± 0.04 M and infers a mass loss of 0.21 ± 0.03 
M between the central RGB and eAGB (the difference in 
initial mass between the RGB and AGB stars is only a few 
xlO -3 M ). The eAGB mass is commensurate with both 
the 0.61 and 0.63 M derived using the Victoria- Regina and 
Dartmouth horizontal branch tracks, respectively (M+09jf|, 
and within the spread of masses (±0.05 M ) which these 
models predict. 

The aforementioned isochrone fits provide the average 
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RGB star we target with another 20 Myr on the RGB. As- 
suming mass loss on the HB itself is negligible, this trans- 
lates to an average mass-loss rate of 10 -8 Mq yr" 1 over 
the remai nder of the RG B. In concurrence with previous 
relations |Catelanl l2009al . and references therein), this is 
rather higher t han the average Reimers' Law mass-loss rate 
|Reimersill975h of 1.9 x 10" 9 M Q yr -1 (assuming Reimers' 
77 = 0.5). It is towards the higher end of the observed range 
of mass- loss rates estimated from chromosphe ric line pro- 
files (e.g. lDupree et al]|2009l ; IVievtes et alj|201ll ). suggesting 
that stronger mass loss at the RGB tip is important. 

One can also use the white dwarf mass to estimate that 
there is only 0.09 ± 0.05 Mq of material for the star to 
lose during the ~9 Myr it spends in its (post-)AGB phases: 
also yielding an average mass- loss rate of 10~ 8 Mq yr -1 . 
Chromospheric mass loss appears con tinuous on the AGB 
and can reach ~6 xlO" 8 M yr" 1 l|Dupree et al.l F2009). 
while the cluster also has several dust-producing stars, 
which typically sustain ~10" 6 M yr -1 for ~10 5 yr (M+09; 
iMcDonald et al.ll201ld V The combination of these processes 
may mean that many stars lose their entire envelope be- 
fore they reach the post-AGB phase, perhaps e ven becom- 
ing A GB-manque stars ('failed' AGB stars; e.g. IO'Connel3 
1999). This would explain the surprisingly-low luminosity 
(L ~ 1500 L Q ) of lj Cen's known post- AGB stars: LEID 
16018 (Fehrenbach's star) and 32029 (VI) jMcDonald et all 
l2011ch . 

Having measured the differential RGB/eAGB mass in 
uj Cen, and confirmed it against existing models, it now be- 
comes possible to repeat the study in other clusters where 
the models are less certain. In particular, this technique 
will be useful in those clusters with higher metallicities and 
higher envelope masses, where the temperature of HB stars 
is less sensitive to stellar mass. Our current observations 
do not allow us to relate metallicity, initial mass, or initial 
abundance to RGB mass loss. For that, we need these ob- 
servations repeated in other clusters. While other clusters 
are less populous than ui Cen, we are limited here by the 
existence of spectra, not by the number of stars. Bespoke 
targetting of stars (particularly on the eAGB), and a se- 
lection of spectral range with more gravity-sensitive lines, 
would allow this study to be repeated to higher accuracy, 
even in much-smaller clusters. 



7 CONCLUSIONS 

We have spectroscopically determined surface gravities for 
66 central RGB and 21 eAGB stars in u Cen, and combined 
these with photometric temperature and luminosity mea- 
surements to calculate the difference in mass between the 
two populations. We find that 26 ± 4% of their mass is lost 
between these evolutionary phases, corresponding to 0.21 ± 
0.03 M© for an initial mass of 0.83 Mq. By implication, this 
limits the mass lost on the AGB to some 0.09 ± 0.05 Mq, 
which may lead to early termination of the AGB and forma- 
tion of AGB manque stars. Our derived HB masses of 0.62 ± 
0.04 Mq compares very well with HB models, which predict 
HB masses of 0.61 - 0.63 M for appropriate (Y w 0.24, 
[a/Fe] » +0.3, [Fe/H] « -1.62) models. This method has 
the potential to provide physical constraints on currently- 
uncertain regimes in modelling horizontal branch stars. 
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